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We report a study of the processes e + e~ — > 777 and e + e~ — > 7/7 at a center-of-mass energy 
of 10.58 GeV, using a 232 fb" 1 data sample collected with the BABAR detector at the PEP-II 
collider at SLAC. We observe 20I5 VI an d 50ly 7/7 events over small backgrounds, and measure 
the cross sections a(e + e" -> 717) = 4.5ih^ ± 0.3 fb and a(e + e^ -> 7/7) = 5.4 ± 0.8 ± 0.3 fb. The 
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and q 2 \F v >(q 2 )\ = 0.251 ± 0.019 ± 0.008 GeV, respectively. 
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I. INTRODUCTION 

The cross section for the reaction e + e~ — ► 7* — » P7, 
where P is a pseudoscalar meson, is given, for energies 
large compared with the P mass mp, by 



dcr e + e 



*p 7 



m:) 



7r 2 a 3 



d cos 9* 



|P P ( S )| 2 (l + cos 2 ^), (1) 



where is the e + e~ center-of-mass (cm.) energy, 8* 
is the angle between the outgoing photon and the in- 
coming electron in the e + e~ cm. frame, and a is the fine 
structure constant. The form factor Fp(q 2 ) describes the 
effect of the strong interaction on the 7* — > 7P transi- 
tion as a function of the four-momentum q of the virtual 
photon; here q 2 = s. 

These transition form factors can be calculated using 
perturbativc Quantum Chromo dynamics (QCD) in the 
asymptotic limit, q 2 »m 2 P [1, 2]: 



-q 2 F P (q 2 ) = V2/p(l 



5 a s (q 2 ) 

3 7T 



(2) 



where fp is the pseudoscalar meson decay constant, and 
a s is the strong coupling. The 7r-meson decay constant 
is known from leptonic ir decays to be about 131 MeV. 
The effective 77 and rf decay constants depend on the 
mixing between the two states, which must be calculated 
from other data [3-8]; for example, the scheme in Ref. 3 
gives /, ; ps f n and f,y ps 1.6/ w [9]. At lower q 2 , how- 
ever, the form factor can only be estimated phenomcno- 
logically. Currently, measurements of e + e _ — > 777 cover 
only the energy region below yfs = 1.4 GeV [10, 11], 
where decays of p(770), u(782), and </>(1020) dominate. 
There are also measurements from J /if)— * 777 and from 
(/), J/tp,^p(2S) ^ rj 1 ^ decays [12]. Space-like 77M7 tran- 
sition form factors have been measured in two-photon 
reactions 77*^ 77 (/) [13-17] up to q 2 ps 20 GeV 2 . These 
q 2 values are not in the asymptotic region, and mea- 
surements at higher q 2 are needed both to establish the 
asymptotic value and to test phenomenological models. 

In this article we present measurements of the reac- 
tion e + e~^ 77W7 at an average e + e~ cm. energy of 
10.58 GeV, corresponding to q 2 = 112 GeV 2 . We recon- 
struct the 77 in the 7r + 7r~7r° decay mode, and the 77' in 
the tt + tt~t] decay mode, where the intermediate 77 state 
decays to cither 77 or 7r + 7r~7r°. From Eqs.(l) and (2) 
and the fp values given above, we expect cross sections 
of 



a(e + 



* 777) « 2.1 fb 
77 7) ps 5.5 fb, 



(3) 



*Also at Laboratoire de Physique Corpusculairc, Clermont- 
Ferrand, France 

tAlso with Universita di Perugia, Dipartimcnto di Fisica, Perugia, 
Italy 

t Also with Universita dclla Basilicata, Potenza, Italy 



which are much smaller than those of many hadronic pro- 
cesses, so we must consider other sources of such events, 
as well as backgrounds, carefully. About 20% of the 
hadronic events in our data are from decays of the T(AS) 
resonance; its branching fraction into 77^7 has not been 
measured, but can be estimated using the relation 



r(T(45) -» ?? (,) 7) r(T(45) 



) 



T(X(1S) -> 77W7) r(T(15) e+e-)' 



(4) 



From the upper limit on the branching fraction 
B(r(lS)-> 77W7) <2.1(1.6)xl0" 5 at 90% CL [12], we 
obtain S (7(45)^ 77 W 7) < 2.5(1.9) xl0~ 8 and a cross sec- 
tion, cT(e+e-^r(45 , )^7/(')7)< 0.026(0.020) fb, well be- 
low the values expected for the mechanism under study. 
Radiative return, e + e~—t r y ISR e + e~ — > JisrV^'\ m which 
there is a high energy photon j ISR from initial-state 
radiation (ISR) off the initial electron or positron, is 
forbidden in single-photon annihilation of the resulting 
e + e~ pair. Double-photon exchange is estimated to have 
a cross section much smaller than in Eqs.(3) [18]. We 
therefore assume all the true e + e~— > 77W7 events in the 
data are due to the processes under study. 

The radiative processes e + e" — > 7 JSf j7r + 7r~7r° and 
e + e~^7j SH 7r + 7r~?7 produce final states identical to those 
for the signals. However, the 7r + 7r _ 7r° and tt + tt~ti mass 
distributions for these processes do not show peaks at the 
77 or 77' masses, and we include this background in the fits 
to the mass distributions. Other sources of non-peaking 
background, such as higher multiplicity ISR events with 



missing particles and e -1 



hadrons events with a high 



energy ir° faking a hard photon, arc reduced to low levels 
in the selection process. 

Background that peaks in the 7/'* 1 mass region arises 
mainly from the ISR processes e + e~^>^ ISR V^"f ISR rii^"f 1 
where V is a vector meson, such as p, uj, <j>, J/ip or Y. 
If the photon from the vector meson decay has low en- 
ergy in the laboratory frame and is lost, the event cannot 
be distinguished from signal. Additional peaking back- 
ground can arise from e + e~ — > VP — > tj^'tt ^/, with or 
without an ISR photon, where V is a vector meson de- 
caying into 7r°7, 777 or 7/7, and P is a ir°, 77, or 77'. These 
backgrounds are estimated from Monte Carlo (MC) sim- 
ulation and data, and are subtracted from the number of 
observed 77^7 events. 



II. THE BABAR DETECTOR AND DATA 
SAMPLES 

Here we analyze a data sample of 232 fb -1 collected 
with the BABAR detector [19] at the PEP-II facility, where 
9.0 GeV electrons collide with 3.1 GeV positrons at a 
cm. energy of 10.54-10.58 GeV. Charged-particle track- 
ing is provided by the five-layer silicon vertex tracker 
(SVT) and the 40-layer drift chamber (DCH), operating 
in a 1.5 T axial magnetic field. The transverse momen- 
tum resolution is 0.47% at 1 GeV/c. Energies of pho- 
tons and electrons are measured with a CsI(Tl) electro- 
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magnetic calorimeter (EMC) with a resolution of 3% at 
1 GeV. Charged-particle identification is provided by 
ionization measurements in the SVT and DCH, and by 
an internally reflecting ring-imaging Cherenkov detector 
(DIRC). Full detector coverage is available over the polar 
angle range 30° < 6* <140° in the cm. frame. 

We simulate the signal processes using a MC genera- 
tor based on Eq.(l). The simulation of ISR background 
processes uses two methods: the Bonncau-Martin for- 
mula [20] for e + er~^ r ) ISR V ', with V — p,uj, (j), J/ip—*rf''l, 
and e + e~ — * ^ ISR u>rf'\ with w—t 7r°7; and the more ac- 
curate approach developed in Ref. 21 where the hadron 
angular distributions are important, for e + e~—t j Isa 3Tr, 
e + e~ — > 7.fSR7r + 7r - f7 and e + e~ — > 7 7SH 4-7r. Since the 
polar angle distribution of the ISR photon peaks near 
0° and 180°, we generate events only over the range 
20° < 6»;<160°, except for e + e~->-/ ISIi r(nS), where the 
ISR photon is generated over the full polar angle range. 
We also simulate non-ISR events of the type e + e _ — *wrf'' 
with to— ► 7r°7. We simulate extra soft-photon radiation 
from the initial state in all cases using the structure func- 
tion method of Ref. 22, with the extra photon energy re- 
stricted such that the invariant mass of the hadronic (plus 
ISR photon) system must exceed 8 GeV/c 2 for non-ISR 
(ISR) processes. We study backgrounds from e + e~—fqq 
using the JETSET [23] package. 

We simulate the detector response, including interac- 
tions of the generated particles with the detector mate- 
rial, using the GEANT4 [24] package, taking into account 
the variation of the detector operating conditions with 
time. We simulate the beam-induced background, which 
may lead to the appearance of extra photons and tracks 
in the events of interest, by overlaying the raw data from 
a random trigger event on each generated event. 



III. EVENT SELECTION 

The initial selection of events requires the presence of 
a high-energy photon with momentum roughly opposite 
to the vector sum of the good-quality charged tracks and 
other photons. The hard photon must have energy in the 
cm. frame E* > 3 GeV; charged tracks must extrapolate 
to the interaction region, have a momentum transverse 
to the beam direction above 100 MeV/c, and have a polar 
angle in the laboratory frame in the region 23° < 9 < 140° 
(38° <9* < 154° in the cm. frame). 

We study the e + e _ — > 777 and e + e~— > 7/7 reactions 
in the Tt + Tr~3j and ir + n~ ir + n~ 37 final states, i.e. we 
use the r/^7r + 7r~7r° decay mode for the former and the 
7/— > 7/7r + 7r~ mode, with 77 — > 77 and rj — ► 7r + 7r~7r , for 
the latter. Since a significant fraction of the events con- 
tain beam-generated spurious tracks and photon candi- 
dates, we select events with at least two (four) tracks 
and at least three photons with energies above 100 MeV 
(50 MeV) for the 2tt37 (4tt37) final state. 

We assume the photon with the highest E* is the recoil 
photon, and consider only the set of two or four tracks 



with zero total charge that has the smallest sum of dis- 
tances from the interaction point in the azimuthal plane. 
We fit a vertex to this set of tracks, which is used as 
the point of origin to calculate all photon angles. We ac- 
cept pairs of other photons as ir° or 77 candidates if their 
invariant mass is in the range 0.07-0.20 GeV/c 2 or 0.45- 
0.65 GeV/c 2 , respectively. For each such candidate, we 
perform a kinematic fit to the selected tracks and pho- 
tons that imposes energy and momentum conservation 
and constrains the candidate ir° or 77 invariant mass. We 
use the x 2 of the kinematic fit {xl^o 1 , xL vl or Xw> 7 ) 
to discriminate signal from background. The simulation 
does not reproduce the shape of the photon energy reso- 
lution function, especially at high energy. Since this dis- 
torts the x 2 distributions, only the measured direction of 
the recoil photon is used in the fit; its energy is a free 
parameter. For events with more than one ir° and/or 77 
candidate, the one giving the lowest x 2 is retained. The 
distribution of X 2 7r7r o 1 for simulated e + e~^777 events is 
shown in Fig. la. There are four effective degrees of free- 
dom, and the distribution shows a long tail due to higher- 
order photon radiation. We also perform the kinematic 
fit without the mass constraint, calculate X2tt3 7 i an d use 
the x 2 difference - x\ v ^ or xLr )7 ~xL 37 ) as a 

measure of the 7r° or r\ reconstruction quality. 

To suppress backgrounds in the e + e~ — > 777 sample 
from events containing kaons and events from multi- 
particle ISR, QED, and e + e~ — > qq processes, while 
maintaining high signal efficiency, we consider events 
with exactly one pair of selected tracks and no more 
than one additional track. Considering the selected 
pair, we require that: i) neither track is identified as 
a kaon; ii) [E\/p\) + (E2/P2) < 1-5, where Ei is the 
EMC energy deposition associated with the i th track and 
Pi is its measured momentum; iii) xj^a-y- xLa 7 < 5 ; 
and iv) the invariant mass of the two charged tracks 
hl-i-K < 415 MeV/c 2 . Requirement (ii) suppresses dielec- 
tron events; requirement (iv) only suppresses background 
events with a 7r + 7r _ 7r° mass > 0.6 GeV/c 2 , but it fa- 
cilitates the extrapolation of the background under the 77 
peak. 

We show scatter plots of xtw^-y versus for the 
selected candidates in the data and the e + e~-^777 sig- 
nal simulation in Figs, lb and lc. A cluster of data 
events is evident near the 77 mass at small values of 
X27r7r° 7 - Figure 2 shows the distribution for data 

events with y 2 . D <20. In order to determine the num- 
ber of events containing a true 77 we perform a binned 
maximum likelihood fit to the spectrum over the 

range 450-650 MeV/c 2 with a sum of signal and back- 
ground distributions. We describe the signal by a sum of 
three Gaussian functions with parameters obtained from 
the simulation, convolved with an additional Gaussian 
smearing function of width ctq = 1.31°'q MeV/c 2 deter- 
mined from high-statistics u;— > ir + ir~ ir° data (see Sec. V). 
The background is a second order polynomial. The line 
on Fig. 2 represents the result of the fit. The fitted num- 
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FIG. 1: Distribution (left) of X L for simulated e + e — > 777^ 7r + 7r 7r°7 signal events. Scatter plots of X2tttt°-i versus the 
7T + 7r~7r° invariant mass for the selected events in data (center) and signal simulation (right). 
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FIG. 2: The 3tt invariant mass spectrum for the e + e~^r]j 
candidates in the data. The dashed curve represents the result 
of the fit described in the text. 



fit to this spectrum similar to that used for the rj signal, 
but over the range 900-1000 MeV/c 2 . The line on Fig. 4 
represents the result of the fit and the fitted number of 
events is iV v =38.1±|-J ± 1.0. 

For the e + e~— >7/'7 reaction in the 47T7r°7 final state 
we require xlmr -/ < 25 and that none of the four charged 
tracks is identified as a kaon. We then search for events 
in which three of the pions are consistent with an rj decay. 
Figure 5 shows the distribution of the 7r + 7r~7r° invariant 
mass (4 combinations per event) for selected candidates 
in the data and e + e~ — > 7/7 signal simulation, with the 
additional requirement that < 1 GcV/ c 2 . Peaks at 
the r\ mass are evident over a modest combinatorial back- 
ground. We select events with at least one combination 
in the range 0.535< M 37r <0.56 GeV/c 2 ; no event in the 
data or simulation has more than one. We fit the 5-7T in- 
variant mass spectrum for the selected data events as for 
the other modes, over the range 900-1000 MeV/c 2 , and 
show the distribution and fit result in Fig. 6. The fitted 
number of events containing a true 77' is 12.0l^±0.3. 



ber of events is N v = 22.7^4 ^ ± 0.6, where the first er- 
ror is statistical and the second is the systematic arising 
from the uncertainty on <jq, variation of the background 
parameters, and using a first or third order polynomial 
background. 

For the e + e~^?/7 reaction in the 7r + 7r _ 37 final state, 
we apply the criteria (i)-(ii) above, the analog of (hi) 
X27r»)7 — X2tt37 < ^' anc ^ a su S ut ly different requirement 
on the two-track mass of M 2lr < 410 MeV/c 2 , the kine- 
matic limit for the two pions from an rf— >7t + 7t~?7 decay. 
Figure 3 shows scatter plots of xItt^ versus M^^n f° r 
selected candidates in the data and the e + e~^rj'j sig- 
nal simulation; a cluster of data events is evident near 
the rj' mass at small values of X^irm- We show the 
spectrum for data events with xi^?^ < 20 in Fig. 4, and 
determine the number of events containing an 7/ with a 



IV. BACKGROUND 

A. e + e~ — ► T77 

We consider both non-peaking and peaking back- 
grounds, where the latter arise from other processes pro- 
ducing true 77 mesons or other mesons whose decays re- 
flect or feed down into the 77 mass region. Figure 2 shows 
that the non-peaking background is small in the 77 mass 
region, but increases sharply toward the upper edge of 
the plot. This is due primarily to the low-mass tail of the 
oj resonance in the ISR processes e + e~ — ► 7 JSfi 7r + 7r~7r , 
and e + e~— >7 /SH 7r + 7r _ 7r 7r°. Our simulation of these pro- 
cesses is tuned to existing data [12], and predicts an M 3w 
spectrum consistent with our selected data both inside 
(excluding the 77 peak) and outside the range of Fig. 2. 
The simulated contributions of other ISR and e + e~^qq 
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FIG. 3: Scatter plots of Xi-n^ vs - Ma,^ for the selected events 
in the data (top) and e + e~ — » ^'7 
signal simulation (bottom). 
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FIG. 5: Distributions of 7r + 7r - 7r invariant mass (4 combi- 
nations per event) for selected events in the data (top) and 
e + e _ ^77'7— >7r + 7r _ 777^7r + 7r _ 7r + 7r _ 7r°7 simulation (bottom). 
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4: The n ir 77 invariant mass distribution for the 
->7/7 candidates in the data with 77'— > TT + ir~ri, 77— ^77. 
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The n n~r] invariant mass distribution for the 
^'7 candidates in the data with 77' — > 7r + 7r - 77, 
'it . The curve represents the result of the fit de- 
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FIG. 7: Scatter plot of xItt77 versus M-^ for the selected data 
events containing an additional photon. 
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FIG. 8: The 3-7T invariant mass spectrum for events in the 
data with xtiw < 25. The solid vertical lines bound the r\ 
signal region; the sideband regions are between these and the 
dashed lines. 



processes to the non-peaking background arc negligible. 

The primary source of peaking background is the set of 
ISR processes e + e~^7 7SR 777, where the 777 comes from 
a p, lo, 4>, or J/ip decay, all of which have been measured. 
We calculate the number of background events using a 
simulation based on the vector meson dominance model 
that includes p, w, and 4> amplitudes with PDG resonance 
parameters [12] and phases of 0°, 0°, and 180°, respec- 
tively, and describes the existing data on the e + e~^rj"f 
reaction in the p-u-(j) mass region [10, 11]. The model 
also includes J/tp production, and predicts a total peak- 
ing background of 2.6±0.5 events. 

The simulation does not include other contributions 
such as decays of excited p, u>, or 4> states, as they are 
unmeasured and expected to be small. As a check, we 
select e + e~— >7 JSR ?77 events explicitly from our data, by 
subjecting any event with an additional photon to a kine- 
matic fit to the 37177 hypothesis. Figure 7 shows a scatter 



TABLE I: The number of selected e + e~ — » jisnrrf events in 
the data in several ranges of the 777 invariant mass compared 
with expectations from the simulation. The first error on each 
expected number is statistical, the second systematic. 



M„ 7 (GeV/c") 


Ndata 


N expect 


0.55-0.95 
0.95-1.05 
1.05-3.05 
3.05-3.15 
3.15-6.50 


25±9 
200±15 
18±12 
31±6 
0.0±1.4 


43±3 ±4 
192±5 ±4 
5±1 ±6 
21±1 ±2 
1±0.4±2 


0.55-6.50 


274±22 


261±5 ±9 



plot of the x 2 °f this fit (xItt-v-y) versus the 3ir invariant 
mass, and Fig. 8 shows the spectrum for events with 
X37r 77 < 25; a strong 77 signal is present. We estimate the 
number of e + e~— >■ JisrVI events by counting the events 
in the signal region indicated in Fig. 8 and subtracting 
the number in the two sidebands. The resulting number 
of events, 274±22, is consistent with the 261±5±9 ex- 
pected from the simulation, where the systematic error 
in the latter is due to experimental uncertainties on the 
input parameters to the simulation. Repeating this exer- 
cise in several different ranges of the 777 invariant mass, 
we obtain the results listed in Table I; data and simula- 
tion are consistent. 

Other possible sources of peaking background are the 
processes e + e~^ VP^ i]ir ~f, where V denotes a vec- 
tor meson, p, u>, or </>, and P is a ir° or 77. The CLEO 
and BES experiments have measured these cross sec- 
tions at yfq* «3.7 GeV [25, 26]; assuming the l/q 4 de- 
pendence of VP form factors predicted by perturbative 
QCD [27], we estimate the e + e~— >r)ir°j cross section to 
be about 3 fb at our cm. energy. The simulated selec- 
tion efficiency is very low due to the additional 7r°, ap- 
proximately 2xl0~ 4 , so we expect only 0.2 background 
events from this source. The corresponding ISR process 
e + e~— > JisrVP can also contribute, and we estimate a 
cross section of about 13 fb, (for 20° < 6* <160°), based 
on our studies of several ISR final states with VP com- 
ponents [28, 29], including 47r, 37t, 37777, 27T77, 2A"7r°, and 
2AT?7. This cross section is relatively large, one-quarter 
of the e + e~— > JjsrVI cross section, but the selection ef- 
ficiency is less than 2xl0~ 5 , so we expect no more than 
0.1 events from this source. 

The e + e~ — > 77 S H J 7 7r °7 arid e + e~ — * 777r°7 events are 
selected about 100 times more efficiently by the 7777 cri- 
teria than by the 777 criteria, and similar factors apply to 
other types of events containing additional pions and/or 
photons. We can therefore make another estimate of 
their overall contribution from the difference between the 
observed and expected numbers of e + e~ — > 7j S «'77 can- 
didates of 13±24 (Table I). Accounting for the -10% 
uncertainty in relative selection efficiencies, we estimate 
< 0.6 such events in our signal peak at the 90% CL. 

ISR production of an T(1S), T(2S), or T(35) rcso- 
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nance could produce a peaking background if the T de- 
cays to 777, since the ISR photon is rather soft. From the 
upper limit on B(T(lS)^rn) of 2.1xl0~ 5 [12], we esti- 
mate that the number of e + e^^^ ISR T{lS)—nrj^^ events 
in our data does not exceed 100. Using the relation 

T(T(nS) -> 777) _ T(T(nS) -> e+e~) 
r(T(lS) -> 777) ~ r(T(15) -► e+e~) ' 

we obtain corresponding limits for T(2S) and Y(3S) of 
50 and 140, respectively. The selection efficiencies for the 
IS, 25, and 35 processes arc below 0.01%, 0.02%, and 
0.08%, respectively, so the total T background does not 
exceed 0.13 events. 

We search for peaking background in the e + e~ — * qq 
process using the JETSET simulation. From 736 mil- 
lion simulated events (corresponding to about twice our 
integrated luminosity) only two events pass the 777 se- 
lection criteria. Only one of them, a K°T' V final state, 
has a 37r invariant mass close to the 77 mass. Since we 
do not expect JETSET to predict rates for such rare 
events correctly, we select e + e~ — > KKrj, 77—5-77 events 
from our data as a check. We perform a kinematic fit to 
the K + K~-f-f hypothesis on all events with at least one 
charged track identified as a kaon, and select events with 
xIkj-j <10- From the 2±30 events found in the data and 
312±14 expected from the simulation, we conclude that 
JETSET overestimates the yield and that this source of 
background is negligible. 

Taking the estimate of the number of peaking back- 
ground events from e + e~— >■ 7/ Sf j777, and considering the 
upper limits on all the other sources as additional system- 
atic errors, we estimate the total peaking background to 
be 2.6±0.8 events. Subtracting this from the number of 
observed events with a true 77, we obtain the number of 
detected e + e~^?77 events: 

N V1 = 20.11^ ± 1.0. 

B. e + e~ — > T7'7 

We estimate backgrounds in the e + e~ — > 7/7 sam- 
ple using similar procedures. The non-peaking back- 
ground is very small for both the 77 — > 77 (see Fig. 4) 
and 77 — ► 7T + 7r~7r° (see Fig. 6) modes. According to 
the simulations, it is dominated by the ISR processes 
e + e~ — ► 7 JSJ? 7T + 7r~7r and e + e~ — * 7 /SH 7T + 7r~7r 7r°. As 
for the e + e~— > 777 process, the simulated non-peaking 
background mass distributions are consistent with those 
observed in data. 

The largest source of peaking background in the sim- 
ulations is the ISR process e + e~ — > 7/ S i?V7, where the 
?7'7 comes mainly from <fi and J/ip decays. These have 
been measured with about 10% accuracy, and we use a 
vector-dominance based simulation similar to that for the 
777 analysis to estimate their contribution. In addition to 
the <j) and J /if), we include contributions from the high- 
mass tails of the p and uj with couplings determined from 



TABLE II: The number of selected e + e~^>'y I sRr]''y events in 
the data in several ranges of the 7/7 invariant mass compared 
with expectations from simulation. The first error on each 
expected number is statistical, the second systematic. 



A/„, 7 (GeV/V) 


-^V data 


N expect 


<1.5 


-2±12 


1.7±0.3±0.2 


1.5-2.0 


6±4 


1.0±0.2±1.0 


2.0-3.0 


2±3 


1.3±0.3±1.3 


3.0-3.2 


97±10 


102±2 ±10 


> 3.2 


3±3 


1.1±0.2±1.1 


Total 


110±17 


107±2 ±10 
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FIG. 9: Distributions of the ir + -K~r] invariant mass for se- 
lected e + e~^ 'yisB.J/ip— > Jisr^'j events in the data (points 
with error bars) and simulation (histogram). 



the measured 77'— >W7 and pj decay widths. We estimate 
a peaking background from this source of 0.3±0.1 events 
in each of the two 77 decay modes. 

We check this prediction by selecting e + e~— > ^ishH'i 
events using a kinematic fit to the 7r + 7r _ ?777 hypothe- 
sis. Selecting events with a X^-nri-yy < 25, wc count signal 
and sideband events in the 7r + 7r~?7 invariant mass dis- 
tributions to obtain numbers of events from this source 
in a set of 7/7 mass intervals. The results from data 
and simulation listed in Table II are consistent. In the 
J / ip mass region these events are practically free of back- 
ground, and we compare the data and simulated n + ir~7] 
invariant mass distributions for events with an 7/7 mass 
in the range 3.05-3.15 GeV/c 2 in Fig. 9. The RMS 
of the distribution is 3.9±0.3 MeV/c 2 in the data and 
3.80±0.06 MeV/c 2 in the simulation. 

To bound peaking background from e + e~ — > rfn 'y, 
e + e~— >-7j SR ?7'7r 7 and other events containing additional 
pions and/or photons, we consider the difference between 
the observed and expected numbers of e + e~ — ► JisrV'i 
candidates in Table II, 3±20. Taking into account the 
factor of 30 difference in selection efficiency, along with 
its 10% systematic uncertainty, we determine that the 
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total background contribution from such processes does 
not exceed 1 event in the rj— > 77 mode or 0.3 events in 
the 77— >7r + 7r~7r° mode. 

From the upper limit B (J (IS)-* 7/7) < 1.6 x 1CT 5 [12], 
we estimate that the number of e + e~— >"f ISR T— ^^isrH'i 
events in our data does not exceed 80, 40, or 110 for the 
IS, 25, or 3S states, respectively. The simulated effi- 
ciencies for such events to pass the e + e~^r/7 selection 
criteria are small, and we estimate that the total T back- 
ground does not exceed 0.03 events for the 77^77 mode, 
and is negligible for the ?7^7r + 7r~7r mode. In the 736 
million e + e~—tqq events simulated by JETSET, we find 
none that passes the 77' 7 selection criteria. 

Considering the upper limits as systematic errors, 
we estimate total peaking backgrounds of 0.3±1.0 and 
0.3±0.3 events in the 27r37 and 47r37 final states, respec- 
tively. Subtracting these from the numbers of observed 
77' events, we obtain a total number of 7/7 events, 

AV7 =49.5+^ ±1.5. 

V. DETECTION EFFICIENCY 

A. e + e~ — > 777 

The detection efficiency determined from the simula- 
tion is Smc = (2.01 ± 0.06)%, where the error includes 
a statistical error and the uncertainty in the value of 
#(?7^7r + 7r~7r ). This efficiency must be corrected to ac- 
count for deficiencies in the simulated detector response. 
We take advantage of the relatively large cross section 
for the ISR process e + e~ — > 7 7SR u;(782)— > j ISIt TT + n~iT , 
which can be selected with very low background [28]. 
The M 37r spectrum for this process is described by 

™= n „ {M )^R c{ M), (5) 

where as w (M) is the Born cross section for e + e~ — » 37r, 
di/dM is the so-called ISR differential luminosity, e(M) 
is the detection efficiency as a function of mass, and R is a 
radiative correction factor (see Ref. 28 for a more detailed 
discussion). The e + e~^37r Born cross section near the us 
mass can be described by a Breit-Wigner function with 
well measured parameters [12]. We calculate the ISR 
luminosity from the total integrated luminosity L and 
the theoretical ISR photon radiator function [30]. The 
radiative correction factor is known with a theoretical 
uncertainty below 1% [22]. We can therefore fit the 37r 
invariant mass spectrum for events passing the criteria 
for this analysis in the us mass region to determine the 
efficiency directly from the data. 

Figure 10 shows this distribution after subtraction of 
the ^0.5% background, estimated from simulation as de- 
scribed in Ref. 28. The fitting function is given by Eq.(5) 
convolved with the simulated detector resolution func- 
tion. There are three free parameters: the efficiency cor- 
rection factor 8 U (e = S u £mc)i the us mass; and ac, an 
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FIG. 10: The tt + tt~iy invariant mass spectrum for data 
events in the us mass region. The curve is the result of the fit 
described in the text. 

ad-hoc Gaussian smearing to account for any resolution 
difference between data and simulation. The curve in 
Fig. 10 represents the result of the fit, which returns: 

8 U = 0.933 ±0.009 ± 0.026, 

a G = 1.3l^±0.3MeV/c 2 , (6) 

where the first error is statistical and the second system- 
atic. The fitted mass is shifted from the nominal value 
by 0.5 MeV/c 2 , consistent with expectations from our de- 
tector simulation. The systematic error in the correction 
factor includes contributions from simulation statistics 
(1.2%), uncertainties on the radiative correction (1%), 
background subtraction (0.2%), and the PDG us width 
(1.5%) and peak cross section (1.5%). The systematic 
error in <jg is due to the uncertainty in the us width. 

Before applying this correction to the e + e~— ^777 effi- 
ciency, we must take into account differences in the dis- 
tributions of any kinematic variables on which the effi- 
ciency depends. Of the many variables studied, three 
show large differences, the photon polar angle # 7 , the in- 
variant mass of the two charged pions M^i and the min- 
imum angle between a charged pion and a photon from 
the 7T° decay 0^ 7 ; in Fig. 11 we compare their distribu- 
tions in simulated e + e~^>^ ISR oj events (solid histograms) 
with those in simulated e + e _ — > 777 events (dashed his- 
tograms). The On-f distribution for the e + e~ — ► 7j Sr cj 
data (dots in Fig. 11) is consistent with that for the sim- 
ulation, but significant inconsistencies are visible in the 
7 (x 2 /dof = 38/14) and M 27r (x 2 /dof = 17/11) distribu- 
tions. To estimate shifts in the efficiency correction due 
to the dependence of the efficiency on these variables, we 
calculate 

r ^J2 J5MT^ P n MC &)> * = (7) 

where P^) is the 6> 7 or M 27r distribution for us-f (777) 
events normalized to unit area, and Xi is the center of 
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FIG. 11: Distributions of the photon polar angle (left), the invariant mass of the two charged pions (middle), and the minimum 
angle between a charged pion and a photon from the n° decay (right) for data (points with error bars) and simulated (solid 
lines) e + e~ '— »7/srW events. The simulated background is shown as the (very small) shaded histograms, and the dashed lines 
show the distributions for simulated e + e~— >rj r -f events. 
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FIG. 12: Distributions of the angle between the two photons from the iv — > 77 or 77—^77 decay (left), and the minimum (middle) 
and maximum (right) energy of the decay photons, for data (points with error bars) and simulated (solid lines) e + e~^ r y IS Rio 
events, simulated background events (small shaded histograms), and simulated e + e~~^r)'"/~ >-K + ir~r) events (dashed lines). 



the z th bin. We obtain the values r 6y = 1.011 ± 0.006 
and rj/j, = 0.973 ± 0.016, from which we calculate the 
efficiency correction 8 v = r0 y rM 2 „5 u) =0.918 ± 0.032, and 
the detection efficiency e = 6^smc = (1.85 ±0.09)%. 



B. e + e — > T7'7 

The simulated efficiency for e + e~— > r/'y—* 27r37 events 
is Smc = (2.91 ± 0.13)%. For this final state we can 
again use the efficiency correction determined from the 
e + e~ — > jisRk> events, taking into account differences in 
the relevant kinematic variables. Considering the same 
set of variables, we find similar results: corrections are 
needed only for 9 1 and with very similar values of 
= 1.016 ± 0.008 and r A / 2x = 0.976 ± 0.013. In ad- 
dition, there arc photon distributions that are different 
for 7T° and rj decays. We show distributions of the an- 



gle between the two decay photons # 77 , and the mini- 
mum and maximum photon energies £^ 7 ,min and £7, max 
in Fig. 12. A disagreement between data and simu- 
lation for e + e _ — > "f IS R^ events is seen in the -E 7jmax 
spectrum (x 2 /dof = 18/11) and we calculate J\e ax = 
1.035±0.016. This correction is not needed for e + e~^?77 
events since their -E 7jma x distribution is very close to that 
for e + e"^7, SR w events. We calculate an efficiency cor- 
rection of 5 V ' = rg rM 27r r E- l max <L> = 0.957 ± 0.037, and a 
detection efficiency of e = 5 T] '£MC = (2-79 ± 0.16)%. 

The simulated efficiency for e + e~^ 7/7— *• 47r37 events 
is (1.05 ± 0.07)%. We estimate an efficiency correc- 
tion for the two additional pions using the ISR process 
e + e~—>"f ISR pr)—>-7r + Tr~rff. We select events in both the 
77— >77 and ?7-^7r + 7r~7r° decay modes with criteria similar 
to those used for the signal. The ir + ir~r] invariant mass 
must be in the range 1.4-1.7 GeV/c 2 where the prj mass 
spectrum is at a maximum, and the invariant mass of 
the p candidate must be in the range 0.64-0.90 GeV/c 2 . 
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From the numbers of selected data and simulated events 
in the two 77 decay modes we determine the double 
ratio 85^ = {N^/N2 1 )data/{N^- K /N2^)MC = 0.98±0.06, 
and we calculate a fully corrected detection efficiency of 
e = 5 57T 5,f£j\fc= (0.99±0.10)%. The ratio of the numbers 
of events selected in the two decay modes, 0.31±0.11, is 
consistent with the ratio of simulated detection efficien- 
cies 0.35±0.04. The total detection efficiency for the two 
modes is (3.78 ± 0.19)%. 



VI. CROSS SECTIONS AND FORM FACTORS 

For each of the two signal processes, we calculate the 
cross section as 



a(e + e 



Pi) 



N Pl 
eL 



R, 



(8) 



where Np 7 is the number of signal events from Sec. IV, e 
is the detection efficiency from Sec. V, L=232 fb -1 is the 
integrated luminosity, and R is a radiative correction fac- 
tor. We calculate R as the ratio of the Born cross section 
for e + e~s-Pj to the total cross section including higher- 
order radiative corrections calculated with the structure 
function method [22]. The simulation requires the invari- 
ant mass of the P7 system Mp 7 >8 GeV/c 2 , for which 
we calculate R = 0.956. The detection efficiency used in 
Eq.(8) is for simulated events with this requirement. The 
value of R depends on the energy dependence of the cross 
section. We use a oc l/q A (sec Eqs. (1) and (2)), and 
investigate the model dependence by recalculating R/e 
under the 1/q 3 and 1/q 5 hypotheses. The relative varia- 
tion is less than 10~ 3 , which we neglect. The theoretical 
uncertainty on R obtained with the structure function 
method does not exceed 1%. We obtain 

a(e + e- -► 777) = 4.5^ ± 0.3 fb, (9) 
a(e + e- -> 777) = 5.4 ±0.8 ±0.3 fb, (10) 

where the first error is statistical and the second system- 
atic. The systematic error is the sum in quadrature of 
contributions from detection efficiency, background sub- 
traction, fitting procedure, and radiative correction. 

The value of R we use does not take into account vac- 
uum polarization, and its contribution is included in the 
results (9)-(10). For comparison with theoretical predic- 
tions, we calculate the so-called "undressed" cross section 
by applying a 7.5±0.2% correction for vacuum polariza- 
tion at 10.58 GeV/c 2 [31], obtaining 



cr(e + e 
er(e + e~ 



* ^T) undressed 
V undressed 



4.21^ ± 0.3 fb, (11) 
5.0tg;? ± 0.3 fb. (12) 



Using Eq.(l) we obtain the values of the 777 and rj'j tran- 
sition form factors at q 2 — 112 GeV 

q 2 \F v {q 2 )\ = 0.229 ±0.030 ±0.008 GeV, (13) 



9 2 |W)I 



0.251 ±0.019 ±0.008 GeV. 



(14) 
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FIG. 13: The magnitudes of the 777 (top) and ^'7 (bottom) 
transition form factors measured in this work (filled circle) 
and by CLEO [13] (triangles) . The shaded boxes indicate the 
ranges of form-factor values calculated according to Eq.(2) 
with the decay constants from Refs. [3-8]. 



VII. SUMMARY 

We have studied the e + e~^?77 and e + e^—nrj'^ pro- 
cesses at an e + e~ cm. energy of 10.58 GeV. We select 



20_ 5 777 and 50_ 7 77' 7 events, measure the cross sections 
and extract the values of the transition form factors at 
9 2 = 112 GeV 2 . 

Since the asymptotic values of the time-like and space- 
like transition form factors are expected to be very close, 
we show our results along with CLEO results for space- 
like momentum transfers [13] in Fig. 13 (we averaged the 
CLEO results obtained in different 77 (77') decay modes). 
The CLEO data rise with increasing q 2 , and are consis- 
tent with the values given by our data points. A precise 
theoretical prediction of the value of the form factor at 
q 2 = 112 GeV 2 is problematic due to uncertainties in the 
effective decay constants, the quark distribution ampli- 
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tudes, and possible gluon content of the rj and 77'. Naively 
taking the decay constants from Refs. [3-8] and calculat- 
ing form factor values according to Eq.(2). we obtain a 
range of values indicated by the shaded boxes in Fig. 13. 
Our data points are at the upper and lower ends of the 
range of predictions for 77 and 77', respectively. The pre- 
dicted ratio of the form factors ranges from 1.6 to 2.3, 
inconsistent with our value of 1.10 ± 0.17. This discrep- 
ancy and the large range of the predictions indicates the 
need for more theoretical input. 
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